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Abstract. We investigate the bulk properties of protoneutron stars in the framework of a relativistic mean 
field theory based on nonextensive statistical mechanics, characterized by power-law quantum distributions. 
We study the relevance of nonextensive statistical effects on the /3-stable equation of state at fixed entropy 
per baryon, in presence and in absence of trapped neutrinos, for nucleonic and hyperonic matter. We show 
that nonextensive statistical effects could play a crucial role in the structure and in the evolution of the 
protoneutron stars also for small deviations from the standard Boltzmann-Gibbs statistics. 
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1 Introduction 

A protoneutron star (PNS) is born immediately after the 
gravitational collapse of a massive star (M « 10-^20 M0). 
During the first seconds of its evolution a PNS is a very 
hot, lepton rich and beta-stable object, with a tempera- 
ture of a few tens of MeV and a lepton concentration typ- 
ical of the pre-supernova matter [1^ . Nearly all of its bind- 
ing energy is in the form of neutrinos trapped inside the 
stellar structure. In particular, neutrinos at the star core 
are basically of the Ve type and have a typically energy 
of about w 200 300 MeV T. Shortly after the PNS 
formation, during the cooling, neutrinos escape through 
the star structure and bring out in few seconds nearly all 
of its binding energy. The total luminosity of this process 
depends from several factors, the most important are the 
total mass of the compact object and the neutrino opa- 
city. During this process, called Kelvin-Helmholtz epoch, 
the PNS evolves in a quasi hydrostatic equilibrium, from 
a lepton rich and hot object to a cold neutrino-free com- 
pact star [3j . From a measure of neutrinos luminosity and 
average energy, which are the most important astrophysi- 
cal observable in the study of the stellar structure, it is in 
principle possible to deduce the total binding energy and, 
therefore, the total baryonic mass of the star [3]. However, 
neutrinos are not only important probes in the study of 
newborn PNS, but have also relevant consequences on the 
chemical composition and in the maximum mass of the 
compact object. The presence of neutrinos in fact alters 
significantly the chemical composition of the star, varying 
the formation threshold of different particle species (in 
general the appearance of non-leptonic negative particles 
is delayed when neutrinos are present). The changes in the 
maximum mass, due to neutrinos trapping, are normally 
greater than those due to finite temperature effects [21131 

HE]. 



From a more theoretical point of view, it appears evi- 
dent that the knowledge of the nuclear EOS plays a cru- 
cial role in the determination of the structure and in the 
evolution of the PNS [B] . The processes related to strong 
interaction should in principle be described by quantum 
chromodynamics. However, in the energy density range 
reached in the compact stars, strongly non-perturbative 
effects in the complex theory of QCD are not negligible. 
The central core of a compact star does not at all resemble 
a quasi-ideal gas of hadrons because strongly dynamical 
correlations are present, including long-range interactions 
[ZlIB] ■ In the absence of a converging method to approach 
QCD at finite density one often turns to (effective and 
phenomenological) model investigations. 

Recently, there is an increasing evidence that the nonex- 
tensive statistical mechanics, characterized by power-law 
quantum equilibrium distributions, can be considered as 
an appropriate basis to deal with physical phenomena 
where strong dynamical correlations, long-range interac- 
tions, anomalous diffusion and microscopic memory ef- 
fects take place pi lTUlfTTlIT^ . A considerable variety of 
physical applications involve a quantitative agreement be- 
tween experimental data and theoretical models based on 
the generalized thermostatistics. In particular, in the last 
years several authors have outlined the relevance of nonex- 
tensive statistical mechanics effects in high energy physics 
and astrophysical problems [ni[Tl[TSl[roi|T71P^|TOl|^|^ 
[^[^[Mll^ . The existence of nonextensive statistical ef- 
fects should strongly affect the finite temperature and den- 
sity nuclear EOS [26,27^. In fact, by varying temperature 
and density, the EOS reflects in terms of the macroscopic 
thermodynamical variables the microscopic interactions in 
different nuclear matter regimes. 

From a phenomenological point of view, in this work 
we plan to investigate the relevance of nonextensive sta- 
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tistical effects on the main physical properties of PNS and 
their related astrophysical implications. At this scope, we 
are going to explore different stages of the PNS evolution. 
The first stage corresponds to an entropy per baryon equal 
to one, in which neutrinos are trapped and strongly influ- 
ence the chemical composition of the PNS. After a short 
time, of about 10 15 s, the temperature of the PNS fast 
rises up until it reaches a value of T ~ 45 80 MeV, it de- 
pends again on the chemical composition [55]. This stage 
is called deleptonization era and corresponds to the maxi- 
mum heating and entropy per baryon {s — 2). This is the 
phase, at high temperature and high baryon density, in 
which the presence of nonextensive effects may alter more 
sensibly the thermodynamical and mechanical proprieties 
of the PNS. Furthermore, let us observe that we are going 
to consider only hadronic composition of the PNS, and we 
do not take in consideration possible formation of a mixed 
hadron-quark phase region or a quark core inside it. This 
choice is principally due to the fact that the presence of 
trapped neutrinos delays the on-set of strange particles 
and also possible formation of quark matter [2j|3]. As a 
consequence, the appearance of a mixed phase or a quark 
core could be shifted at the end, or nearly at the end of 
the deleptonization era |29) . 

The paper is organized as follows. In Section ^ we 
briefly introduce the basic formalism of the nonexten- 
sive statistics. In Section [3l we study the nonextensive 
hadronic EOS in the context of a relativistic mean field 
model, in presence and in absence of neutrinos and hyper- 
ons. In Sectional we show the main results in the particles 
concentrations and bulk stellar properties in the presence 
of nonextensive statistical effects. Finally, we summarize 
our conclusions in Section O 



2 Basic assumptions in nonextensive statistics 

Starting point of the nonextensive statistical mechanics, 
introduced by Tsallis, is the following definition of a q- 
deformed entropy functional [9.,10..11j: 



9-1 



1 



/(x)dr2-i (1) 



where /(x) stands for a normalized probability distribu- 
tion, X and df2 denoting, respectively, a generic point 
and the volume element in the corresponding phase space 
(here and in the following we set the Boltzmann and the 
Planck constant equal to unity). The generalized entropy 
has the usual properties of positivity, equiprobability, con- 
cavity and irreversibility, preserves the whole mathemat- 
ical structure of thermodynamics (Legendre transforma- 
tions) . The real parameter q determines the degree of non- 
additivity exhibited by the entropy form ([T]) and in the 
limit g — > 1, becomes additive and reduces to the stan- 
dard Boltzmann- Gibbs entropy. 

A peculiarity of the Tsallis generalized thermostati- 
stics is that if we have two statistically independent sub- 
systems Ae B, described, respectively, by individual proba- 
bility density /(^^ and /(^^ and we call /(^+^)(x^, x^) = 



/'-'^•'(xyi) /*^^^(xb) the joint probability density of a com- 
posite system A + B, the nonadditive (nonextensive) char- 
acter of Sq is summarized in the relation [llj 

^,[/(^+^']= ^,[/(^'] + ^J/(^']) 

+ (l-g)5,[/(^)]5J/(^)]. (2) 

In the limit g — > 1 , the above equation reduces to the well- 
known additivity (extensivity) relation of the Boltzmann- 
Gibbs logarithmic entropy. Here, the word nonextensive 
should be associated with the fact that the total energy 
of long-range- interacting mechanical systems is nonexten- 
sive, in contrast with the case of short-range-interacting 
systems, whose total energy is extensive in the thermody- 
namical sense |11|. 

Second crucial assumption on nonextensive statistics 
is the introduction of the g-mean value (or escort mean 
value) of a physical observable ^(x) [51 [TU1[TT| : 

/A(x) [/(x)]^rfJ? 



(3) 



The probability distribution can be obtained maximiz- 
ing the measure Sq under appropriate constraints related 
to the previous q-mean value definition. In this context, 
it is important to observe that the Tsallis classical distri- 
bution can be seen as a superposition of Boltzmann dis- 
tributions with different temperatures which have a mean 
value corresponding to the temperature appearing in the 
Tsallis distribution. The nonextensive q parameter is re- 
lated to the fluctuation in the temperature and describes 
the spread around the average value of the Boltzmann 
temperature |16j . 

Following the above prescriptions, it is possible to ob- 
tain the associate quantum mean occupation number of 
particles species i in a grand canonical ensemble. For a 
dilute gas of particles and for small deviations from the 
standard statistics {q ~ 1), it can be written as 30, 3l] 

where /3 = 1/T and the sign (±1) is for fermions and 
bosons respectively. Furthermore, in Eq.Q, following Ref. 
[?T] . for q > 1, we have 



,(^) 



l[l 

whereas, for q < 1, 

[1 
[1 



{q - l)x 
(1 - q)x_ 



1/(1-9) 



if a; > 
if a; < 



(5) 



,(^) 



]i/(?-i) 

11/(1-9) 



if o; < 
if a; > 



(6) 



(9-1) 
(1 - <1>] 

Naturally, for q ^ 1 the above quantum distribution re- 
duces to the standard Fermi-Dirac and Bose-Einstein dis- 
tribution. Let us observe nonextensive statistical effects 
vanishes approaching to zero temperature. On the other 
hand, the nonextensive statistics entails a sensible differ- 
ence on the power-law particle distribution shape in the 
high energy region with respect to the standard statistics. 
Hence, nonextensive effects could play an important role 
in the finite temperature and high baryon density PNS 
evolution. 
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3 Nonextensive hadronic equation of state 



example, for g > 1 and /3{E* - \ fi*\) > 0, we have 



In this Section, we study the nonextensive hadronic EOS 
in the framework of a relativistic mean field theory in 
which baryons interact through the nuclear force medi- 
ated by the exchange of virtual isoscalar-scalar (a), isosca- 
lar-vector (w) and isovector-vector (p) mesons fields (32l 
[551134) . In our analysis we include all the baryon octet 
in order to reproduce the chemical composition of the 
PNS at high baryon chemical potential. We also take into 
account of leptons particle by fixing the lepton fraction 
Yl = Ye + Y^^ = {pe + Pi,J/pB, where pe, Pi,^ and ps 
are the electron, neutrino and baryon number densities, 
respectively. This is because, in the first stage of PNS 
evolution, electrons and neutrinos are trapped inside the 
stellar matter and, therefore, the lepton number must be 
conserved until neutrinos escape out of the PNS 

The Lagrangian density can be written in term of the 
hadronic [34[|55] plus leptonic component, as follow: 

jC-tot ^ Ch + Ci -ipBiijf^d^ - {Mb ~ g^Bcr) 



-gujBip.'^^ - gpBi'^T 



(7) 



where the sums over B and I are over the baryon octet and 
lepton particles, respectively. The field strength tensors 
for the vector mesons are given by the usual expressions 
= df,uj^ - dyUJf_,, Gf,^ = df,py ~ d^p^, and U{a) is a 
nonlinear potential of a meson 



U{a) = laa^ + \ba^ 
^ ^ 3 4 



(8) 



usually introduced to achieve a reasonable compression 
modulus for equilibrium nuclear matter. 

The field equations in a mean field approximation are 

(^7^9'^ - (Af - g^B(T) - g..B7°oj - gpBl'^T^p)^ = , (9) 
m^cr + acr^ + 6cr^ = g„B < V'V' >= g^BPS , (10) 
mluj = gi^B < -ipj^ip >= guBPB , (11) 



"f^pP = gpB < -07 TaVj >= gpBPi : 



(12) 



where a = (a), uj — (w") and p — are the nonvan- 
ishing expectation values of mesons fields, pj is the total 
isospin density, pB and ps are the baryon density and the 
baryon scalar density, respectively. They are given by 



PS 



i=B 



d^fc M* 



(13) 



(27r)3 E* 



[n!{k) + n^Hk)], (14) 



where ni{k) and ni{k) are the q-deformed fermion parti- 
cle and antiparticle distributions given in Eq.slH])-®; for 



ni{k) 



1 



(15) 



[l + {q-l)f3{E*{k)-fl*)]m,-l) + l 

"'^^^ = [l + (g-l)/?(i?*(fc) + M*)]i/('-i) + l ■ ^^^^ 
The nucleon effective energy is defined as Ei*{k) — 



k^ + Mi*, where Mi* = M, —gaBcr. The effective chem- 
ical potentials /x* are given in terms of the meson fields as 
follows 



Pt 



fj-i - g^B^^ - TsiBgpBP , 



(17) 



where fii are the thermodynamical chemical potentials 
{fii = de/dpi). 

The further conditions, required for the /3-stable chem- 
ical equilibrium and charge neutrality, can be written as 

PA = Pe« = Peo = Pn , (18) 
Ps- = Pe- = Pn + Pe , (19) 

Pp = PE+ = Pn- Pe ; (20) 

Pp + Pe+ - Pe- - Pe- - Pe = . (21) 

In the case of trapped neutrinos, the new equalities are ob- 
tained by the replacement of pe Pe — Pu^ ■ The total en- 
tropy per baryon is calculated using s — [Sb + Si) / {T ps), 
where Sb ^ Pb + ^b - J2i=B PiPi = Pi + - 

PiPiy s-nd the sums are extended over all the baryons 
and leptons species. 

The thermodynamical quantities can be obtained from 
the thermodynamic potential in the standard way. More 
explicitly, the baryon pressure Pb and the energy density 
eB can be written as 



Pf 



(27r)3i?*(fc) 



£B 



-U{a) + \ml,? + , (22) 

2E / 7KZ^E:ik)[nlik)+n1{k)] + -m^ 



(27r)3' 
1 

2' 



(23) 



The numerical evaluation of the above thermodynami- 
cal quantities can be performed if the meson-nucleon and 
meson-hyperon coupling constants are known. Concerning 
the meson-nucleon coupling constants (ga-N, g^N, gpN), 
they are determined to reproduce properties of equilib- 
rium nuclear matter such as the saturation densities, the 
binding energy, the symmetric energy coefficient, the com- 
pression modulus, and the effective Dirac mass at satura- 
tion. Because of a valuable range of uncertainty in the 
empirical values that must be fitted, especially for the 
compression modulus and for the effective Dirac mass, 
in the literature there are different sets of coupling con- 
stants. Here and in the following, we focus our investiga- 
tion by considering the so-called GM3 parameter set (even 
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if comparable results can be obtained in other parameter 
sets) [341135] . The implementation of hyperon degrees of 
freedom comes from determination of the corresponding 
meson-hyperon coupling constants that have been fitted 
to hypernuclear properties and their specific values are 
taken from Ref.[3S] for the GM3 parameter set. 

4 Results and discussion 

As briefly mentioned in the Introduction, our analysis fo- 
cus on the relevance of possible nonextensive statistical 
effects during the first PNS evolution phases. We can ide- 
ally divide the evolution into three phases. The first, at 
the beginning, in which neutrinos are trapped and the en- 
tropy per baryon is assumed fixed to s = 1 and = 0.4. 
A second phase, after about 10-^15 sec, which corresponds 
to the maximum heating of the star and neutrinos are free 
(s = 2,Y,^^ = 0). Finally, a third phase of cold-catalyzed 
PNS (s = 0,y^^ = 0) [11133. Regarding the relevance of 
nonextensive statistical effects, the most important phase 
corresponds to the maximum heating, in which the pres- 
ence of nonextensive statistical effects may play a crucial 
role in the determination of the PNS chemical composition 
and related thermodynamical proprieties. In this work we 
limit ourselves to consider only a small variations from the 
standard statistics (from q = 0.97 to q = 1.03). 

In Fig. [TJ we show the temperature as a function of 
the baryon density (in units of the saturation nuclear 
density po = 0.153 fm^"^) and for different values of the 
nonextensive parameter, in absence {np) and in presence 
(npH) of hyperons. We limit our analysis in the first two 
phases: in the upper panel, the first leptonic rich state 
(s = Iji^L = 0.4) and, in the lower panel, the maxi- 
mum heating phase (s — 2,Y^^ — 0). Indeed in the cold- 
catalyzed phase (s = 0,Yi,^ — 0), the temperature is very 
low (fews MeV), and nonextensive statistical effects may 
be neglected. In both previous cases, we observe a reduc- 
tion in temperature in presence of a sub-extensive stati- 
stics {q < I) and a general increase for g > 1. This effect 
is more remarkable when hyperons are present and for 
higher values of entropy for baryon. For example, in the 
maximum heating phase, in presence of hyperons and for 
g = 1, the temperature is about T = 37 MeV at baryon 
density pB = 5po, whereas, for q = 1.03 and q = 0.97, it 
is approximately equal to T = 51 MeV and T = 21 MeV, 
respectively. Note also that, when hyperons are present, 
for s = 1 and Yl = 0.4, the system evolves in a quasi 
isothermal configuration above pB — 2.5-^3po- The dif- 
ferent behavior in the stellar temperature have important 
consequences in the PNS evolution and in its particles 
concentration. Finite temperature properties of matter at 
high density influence the diffusion of neutrinos, being the 
neutrino mean free paths strongly temperature dependent 
[Hl37j. In particular, neutrino opacity is very sensitive to 
the inner temperature (in general proportional to T^) and, 
therefore, this would affect sensibly the cooling of the 
PNS, making it longer when g > 1, and shorter when 
q < I. This matter of fact could have important conse- 
quences on the neutrino luminosity, because its drop is 



T[MeV] 




Fig. 1. Temperature as a function of the baryon density (in 
units of the saturation nuclear density po) for different values 
of q, entropy per baryon and neutrino fraction (upper panel: 
s — 1,Yl = 0.4 and lower panel: s — 2, V^^ = 0). The labels 
np and npH stand for nucleons and nucleons plus hyperons. 



associated with the end of this cooling process [3]. Con- 
sequently, an alteration in the Kelvin-Helmholtz epoch, 
that does not correspond to the predict direct or modified 
URCA process, could be an indication of nonextensive sta- 
tistical effects. 

In Fig. [H is reported the dependence of the pressure 
from the baryon density, for different values of q, in the 
initial phase: s = 1,Yl — 0.4 (upper panel) and in the 
maximum heating phase: s = 2 and Y^^ = (lower panel). 
With the appearance of hyperons, around pB = 3poj we 
have a general softening of the EOS. However, due to the 
low temperature achieved in this phase (see Fig.[T]), nonex- 
tensive statistical effects do not change significantly the 
total pressure of the PNS. The situation is somewhat dif- 
ferent when we analyze the maximum heating phase. In 
such a condition, the temperature is higher and nonexten- 
sive statistical effects are more relevant, especially when 
hyperons are present. 

As it is well known, the softening of the EOS, due to 
the appearance of additional fermionic degrees of freedom 
in the form of hyperons, leads to higher central densities. 
This matter of fact is, however, influenced from the pres- 
ence of nonextensive statistical effects. To better focalized 
this aspect, in Fig. [31 is plotted the central baryon den- 
sity pc in the PNS core corresponding to a total baryon 
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Fig. 2. Total pressure as a function of baryon density (in units 
of the saturation nuclear density) for different values of q in the 
lepton rich case: s — 1 and Yl = 0.4 (upper panel) and in the 
maximum heating phase: s = 2 and Y,^^ = (lower panel). 
The labels np and npH stand for nucleons and nucleons plus 
hyperons matter, respectively. 



mass Mb (in units of the solar mass Mq). The case q > I 
(g < 1) implies a general reduction (increase) of the cen- 
tral density at fixed total baryon mass. This effect is em- 
phasized in the maximum heating condition (lower panel). 
In particular, for g > 1 and until Mb ~ 1.7 Mq, the cen- 
tral densities of hyperons stars (npH) are lower than the 
ones corresponding to the standard (g = 1) nucleons-only 
stars {np). This feature, together the variation of the tem- 
perature, can be very relevant in the determination of the 
specific heat of the stellar matter and, as a consequence, 
on the neutrino diffusion [3]. 

In addition to the above considerations, we investi- 
gate the lepton chemical potentials in different PNS con- 
ditions because they strongly influence the deleptoniza- 
tion process [HIS]. At this scope, in Fig. |31 we show the 
lepton chemical potentials in lepton rich matter (s = 1 
and Yb = 0.4) in absence (np) and in presence {npH) of 
hyperon degrees of freedom (higher curves for the elec- 
tron chemical potential //e and lower curves for the neu- 
trino chemical potential ^.v^)- For clarify, we report only 
the case q > 1 compared to the standard {q = 1) re- 
sults. It is interesting to observe that for a nucleons-only 
EOS, super-extensive statistical effects imply a general re- 
duction of the electron and neutrino chemical potentials 



PclPa 




Mfi[M„] 



Fig. 3. Central baryon density pc (in units of the saturation 
nuclear density) corresponding to a total stellar baryon mass 
Mb in the lepton rich phase (upper panel) and in the maximum 
heating phase (lower panel). The labels np and npH stand for 
nucleons and nucleons plus hyperons matter, respectively. 



respect to the standard case. Otherwise, in presence of 
hyperons, a sensible reduction of the electron chemical 
potential is not accompanied by a reduction of the neu- 
trino chemical potential which retains very similar values 
to the standard case at high baryon density. Such a fea- 
ture could have important consequence on the diffusion of 
the electron neutrinos inside the PNS. 

Moreover, the neutrino mean free paths and the mat- 
ter specific heat depend sensitively on the composition; 
under degenerate conditions even modest changes to the 
composition significantly alter the neutrino scattering and 
absorbtion mean free paths. It is, therefore, relevant to 
investigate how nonextensive statistical mechanics influ- 
ences particle compositions in different PNS epochs. In 
Fig. [5l we report the particle concentrations for s = 1 
and Yl — 0.4 in absence (upper panel) and in presence 
(lower panel) of hyperons for different values of q. It is 
well known that the presence of trapped neutrinos signifi- 
cantly alter the protons and the electrons abundance and 
strongly influence the threshold of hyperons formation [21 
|31[S]. In absence of hyperons, nonextensive statistical ef- 
fects do not play a signiflcantly role, in fact the particle 
concentrations are almost the same. The situation changes 
when we include hyperon degrees of freedom. The presence 
of sub-extensive effects {q <1) slightly lowers the neutri- 
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Fig. 4. Lepton chemical potentials in lepton rich matter as a 
function of the baryon density (in units of the nuclear satura- 
tion density). The labels np and npH stand for nucleons and 
nucleons plus hyperons matter, respectively. 



nos concentration, while increases the neutrons and the 
electrons ones. Moreover, hyperons start later and their 
concentration are in general decreased, excepted for the 
A particles ratio, which becomes greater after ps > 4 po 
respect to the standard case. Otherwise, in case of super- 
extensive statistical effects {q > 1), we can observe a gen- 
eral reduction in the neutrons and the electrons fractions 
and a small increase of the neutrinos and the hyperons 
concentrations, with important effects on the softening of 
the EOS. 

In the Fig. El particle concentrations for the maximum 
heating phase (s = 2, Y^^ = 0) are reported. In this condi- 
tion, due to the higher temperature achieved in the stellar 
matter, nonextensive statistical effects become more rel- 
evant and, consequently, particle concentrations change 
significantly. When q < 1, we observe a reduction in the 
protons and the electrons concentrations and an increase 
of the neutrons fraction. Whereas, in the super-extensive 
case {q > 1), we have a lower neutrons fraction and an in- 
crease of the protons and electrons concentrations. When 
hyperons are included. Fig. [SI lower panel, we have two 
main consequences. Firstly, with the absence of neutrinos, 
the hyperons on-set is shifted at low baryon densities, be- 
low 2 po- Secondly, as a consequence, the electrons and the 
protons concentrations decrease sensibly with respect to 
the initial lepton- rich regime (s = 1, = 0.4). Therefore, 
the absence of neutrinos in the stellar matter implies a 
strong softening of the EOS. In presence of sub-extensive 
statistics {q < 1), hyperons start later with respect to the 
standard case {q — 1), and have in general a bigger con- 
centration at high baryon density. The other way round 
takes place for super-extensive statistics (q > 1). 

To better understand the role of nonextensive statisti- 
cal effects on hyperons formation in the PNS core, in Fig. 
[71 we report the total hyperons concentration (strange- 
ness per baryon) as a function of the stellar baryon mass 
(in units of Mq) in the lepton rich (upper panel) and in 
the maximum heating (lower panel) epoch. Although, for 
q > 1 , hyperons are present at lower baryon densities with 



0.1- n 

0.6^ 



0.5: 




Fig. 5. Particle concentrations Yi without (upper panel) and 
with hyperons (lower panel) as a function of the baryon density 
for s = 1 and Yl = 0.4. 



respect to the standard case {q — 1) and a greater hyper- 
ons concentration at lower baryon masses is present, we 
have a significant reduction at higher baryon masses, es- 
pecially in case of the maximum heating condition. For 
q < 1, we have instead the opposite effect: a reduction of 
the hyperons fraction at low baryon masses and an en- 
hancement at high baryon masses. These features, prin- 
cipally due to the behavior of A particles concentration 
as a function of the baryon density (see Fig.s [Sj and [H]), 
could imply relevant phenomenological consequences on 
the evolution of the PNS. In fact, it is known that hyper- 
ons significantly increase the neutrino scattering and ab- 
sorption cross sections j4^. Furthermore, the central den- 
sities of hyperons stars become progressively larger than 
that of purely nucleon stars and the evolution timescale of 
hyperons stars results to be slightly larger because of the 
smaller mean free path of hyperonic matter. Larger central 
densities and higher electron neutrino energies, reached in 
hyperonic PNS, increase the neutrino opacity, temporar- 
ily reducing the loss of neutrinos from the stellar core and 
allow to sustain a higher luminosity at late times. 

Finally, in the Table [TJ we report the maximum grav- 
itational (baryonic) masses (in units of Mq) and the cor- 
responding values of radius and central baryon density, 
for different PNS configurations. We have taken in consi- 
deration different values of the nonextensive g-parameter. 
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Yi 




Fig. 6. The same of Fig. [S] for the maximum heating phase 
(s = 2 and Y^^ = 0). 



entropy per baryon and leptons concentration, in absence 
(np) and in presence (npH) of hyperons. 

As we have already remarked, nonextensive statisti- 
cal effects result to be most relevant when hyperons are 
present and the maximum heating phase is achieved (s = 
2, =0). In fact, when the entropy per baryon is equal 
to one and the matter is lepton rich, the difference in the 
maximum mass is very small both with and without hy- 
perons. Moreover, immediately after the deleptonization, 
nonextensive statistical effects become more important, 
especially in presence of hyperons. The maximum mass 
is in general decreased when we consider a sub-extensive 
statistics {q < 1), and increased when q > 1. 

Furthermore, in the presence of hyperons, when the 
stellar core contains non-leptonic negative charges, the 
maximum masses of neutrino-trapped stars result to be 
significantly larger than for low temperatures and for lep- 
ton poor matter. Hence, there exists a window of initial 
masses for which the star becomes unstable to gravita- 
tional collapse during deleptonization and a black hole can 
take place [21IHIS]. We can see that the formation of such 
a metastable phase strongly depends on the presence of 
nonextensive statistical effects and the window of metasta- 
bility grows with the value of the nonextensive index q. In 
particular, for g > 1, we still have a very large value of the 
maximum baryonic mass in the maximum heating phase 
(s = 2, Yy_. = 0), significantly larger than the one in the 



Ps/pb 




Fig. 7. Hyperons concentration Ps/pB (strangeness per 
baryon) as a function of the baryon mass Mb in the lepton 
rich (upper panel) and in the maximum heating phase (lower 
panel) . 



cold catalyzed condition (where nonextensive statistical 
effects do not play any role). Therefore, for nucleons-only 
stars, a black hole could form only after the core bounce, 
because the maximum mass supported by neutrino-free 
stars is bigger than that supported for neutrino-rich case. 
Whereas, in presence of hyperons, a black hole can take 
place also after the deleptonization era because of the re- 
alization of a metastable phase, which becomes more rel- 
evant in presence of super-extensive statistical effects. 



5 Conclusions 

We have investigated the physical properties of the PNS 
in the framework of a relativistic mean field theory based 
on nonextensive statistical mechanics, characterized by 
power-law quantum distributions. We have studied the fi- 
nite temperature EOS in /3-stable matter in absence and 
in presence of hyperons and trapped neutrinos. From a 
phenomenological point of view, we have considered the 
nonextensive index g as a free parameter, even if, in prin- 
ciple, it should depend on the physical conditions inside 
the PNS, on the fluctuation of the temperature and be 
related to microscopic quantities (such as, for example, 
the mean interparticle interaction length). In this context, 
let us remember that, in the diffusional approximation, a 
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Table 1. Maximum gravitational (baryonic) masses Afmax (in 
units of M©) and corresponding values of radius R and central 
baryon density pc in absence (np) and in presence (npH) of 
hyperons for different values of the nonextensive parameter 
q. The results arc reported for different values of entropy per 
baryon and leptons concentration. For completeness, the values 
for the cold-catalyzed phase (s = 0, Y^^ = 0) for g = 1 are also 
reported (nonextensive statistical mechanics does not play any 
role in this last rc\e;im(>). 





np 


q 


s^1,Yl= 0.4 


s = 2, y^,^ = 


0.97 


M„,ax = 1.96 (2.20) 
R= 10.52 km 
Pc = 7.33 Po 


Af^ax = 2.06 (2.34) 
R=11.37 km 
pc = 6.50 Po 


1.00 


M^ax = 1.97(2.21) 

R=10.69 km 
pc — 7.23 po 


Af^ax = 2.09 (2.38) 
R=11.80 km 
Pc = 6.17 po 


1.03 


Afniax = 1.98 (2.23) 
R=10.81 km 
Pc = 7.15 Po 


Mn,ax = 2.15(2.47) 
R= 12.39 km 
Pc = 5.75 Po 




s = 0, = 


1.00 


Mn,ax = 2.05 (2.39) 
R=ll.ll km 
Pc = 6.91 po 




npH 


q 


s=1,Yl= 0.4 


s = 2, F,, = 


0.97 


Afniax = 1.75 (1.93) 

R=ll.ll km 
Pc = 6.84 Po 


A/n.ax = 1.55 (1.69) 
R=12.39 km 
Pc = 5.59 po 


1.00 


AT^ax = 1.76 (1.94) 
R=11.41 km 
Pc = 6.62 po 


Afn,ax = 1.59 (1.74) 
R=12.75 km 
Pc = 5.44 Po 


1.03 


Afniax = 1.78 (1.96) 
R=11.56 km 
Pc = 6.51 Po 


Afniax = 1.71 (1.90) 

R=13.34 km 
Pc = 5.05 Po 




s = 0, = 


1.00 


Afn,ax = 1.57(1.76) 
R=12.35 km 
Pc = 5.66 Po 



value q ^ 1 implies the presence of an anomalous dif- 
fusion among the constituent particles (the mean square 
displacement obeys to a power-law behavior < x"^ >oc 
with a 7^ 1). 

We have shown that nonextensive statistical eflFects 
could play a crucial role in the structure and in the evo- 
lution of the PNS also for small deviations from the stan- 
dard Boltzmann-Gibbs statistics. As expected, nonexten- 
sive statistical effects result to be particularly important 
during the maximum heating phase (s = 2,Yi,^ = 0), 
while become less relevant during the initial lepton rich 
state (s = 1,Yl = 0.4) and negligible in the cold cat- 



alyzed phase (s = 0, V^^ = 0), due to the low temperatures 
achieved. 

We have studied the relevance of nonextensive statis- 
tical effects; i) in the temperature behavior as a function 
of the baryon density, ii) in the softening of the EOS and, 
consequently, in the central baryon densities reached in the 
PNS core at fixed baryon mass, iii) in the lepton chemical 
potentials, iv) in the particle concentrations, in the hyper- 
ons formation and in the strangeness per baryon at fixed 
total baryon mass, v) in the maximum gravitational and 
baryonic masses. Such a variation of physical quantities, 
respect to the standard case, can imply important conse- 
quences on the determination of the matter specific heat, 
on the neutrino mean free path inside the stellar core and, 
consequently, in the neutrino opacity and luminosity. 

We have considered both cases of sub-extensive {q < 1) 
and super-extensive (q > 1) statistical effects which entail 
a sensible difference on the power-law particle distribution 
in the high energy region. When the entropic q parameter 
is smaller than one, the energy tail of the particle distri- 
bution is depleted, otherwise, when q is greater than one, 
the energy tail is enhanced. In the PNS context, the phys- 
ical meaning of the difference between a sub-extensive or 
a super-extensive statistical behavior is reflected in differ- 
ent and well distinguishable phenomenological PNS prop- 
erties. For g < 1, we have a reduction in temperature at 
fixed baryon density with respect to the standard case 
(q = 1). Especially in the maximum heating phase, the 
EOS becomes slightly softer and higher central baryon 
densities at fixed total baryon mass are reached, infiuenc- 
ing the neutrino diffusion during the deleptonization pro- 
cess. In the case of sub-extensive effects, hyperons start 
later but have a bigger concentration at high baryon den- 
sity, allowing to sustain a higher neutrino luminosity at 
late times. The other way round takes place in the case 
oi q > 1. We have an increase of the temperature as a 
function of the baryon density and lower central densities 
at fixed baryon mass are rcaciied. The hyperons on-set 
is shifted at lower baryon densities and a greater hyper- 
ons concentration at low baryon masses is present. On the 
other hand, a significant rediiction of the hyperons concen- 
tration at high stellar masses take place, contributing to a 
lower luminosity at late times. Moreover, we have shown 
that, in presence of supcr-cxtcnsivc statistical effects and 
hypcron degrees of freedom, it is favored the realization 
of a metastable phase, with an enhancement of a possible 
black hole formation after the deleptonization era. 
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